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Realistic vacuum tttv interactions are employed to investigate thermal tt + tt~ emission spectra from 
the late stages of heavy-ion reactions at ultrarelativistic energies. Hadronic in-medium effects, in- 
cluding many-body p-meson spectral functions used earlier to describe the dilepton excess observed 
at CERN-SPS energies, are implemented to assess resulting modifications in relation to recent mea- 
surements of it + -k~ invariant-mass spectra by the STAR collaboration in p-p and Au-Au collisions at 
y / 5ivjv=200 GeV. Statistical model estimates for the p / 71 " - and K* /K ratios close to the expected 
thermal freezeout are also given. 



I. INTRODUCTION 

The unambiguous identification of hadron modifica- 
tions in hot and/or dense matter constitutes one of 
the major objectives in modern nuclear physics [1,2]. 
Hadronic states whose vacuum properties are related 
to the spontaneous breaking of chiral symmetry, are of 
particular interest as changes in their spectral distri- 
bution reflect precursor phenomena towards the chiral 
phase transition. A versatile environment for these in- 
vestigations is provided by high-energy heavy-ion colli- 
sions, as different center-of-mass (CM) energies have been 
shown [3] to probe large regions in temperature and den- 
sity of the QCD phase diagram. 

On the one hand, dilepton observables [4-6] have been 
proved to be powerful tools to extract medium modifi- 
cations of vector-mesons [2] during the evolution history 
of the hot and dense fireballs formed in heavy-ion reac- 
tions. However, the relative signal weakness (paramet- 
rically suppressed by a 2 compared to hadronic sources) 
entails appreciable experimental uncertainties, which, as 
of now, do not allow for a unique theoretical interpre- 
tation of the results. On the other hand, hadronic fi- 
nal state measurements of strongly decaying resonances 
have thus far not revealed significant medium effects [7,8] . 
This is believed to be mostly due to the diluteness of the 
matter from which the decay products are able to reach 
the detectors undistorted (unlike in the dilepton case). 
In this context, recent (preliminary) data [9] on tt + tt~ 
invariant-mass distributions from GeV nu- 

clear collisions by the STAR collaboration are of partic- 
ular interest. Assuming Breit-Wigner parametrizations 
of the underlying resonance structures (p, w, fo, etc.), 
a downward shift of ~ 30 MeV for the peak position in 
the 7r + 7r~ distribution from p-meson decays has been re- 
ported [9] when going from p-p to 40-80% central Au-Au 
collisions. In addition, a shift of -40 MeV relative to a 
nominal p-meson mass of 770 MeV has been extracted 
from the p-p data. 

Two recent articles have specifically addressed these 
observations (see also Ref. [10] for an earlier study). In 
Ref. [11], a number of hadronic s- and i-channel ex- 
changes between p-mesons and the surrounding mat- 
ter close to freezeout have been estimated to produce 



a ^50 MeV downward mass shift, mostly originating 
from scalar f-channel exchanges with anti-/baryons (in 
line with " Brown- Rho" scaling [12]). In Ref. [13] a more 
generic analysis of the interplay between phase-space, 
mass-shift and width effects has been carried out, with 
the latter two implemented via schematic selfenergies in 
the underlying p-meson spectral function. 

In the present article we go beyond these analyses in 
several respects. We will assess ir + ir~ spectra based on 
microscopic 7T7t interactions that accurately describe per- 
tinent free scattering phase shifts in S-, P- and D-waves, 
including hadronic medium modifications through finite- 
temperature effects on intermediate two-pion states. For 
the p meson we will also use more advanced in-medium 
many-body spectral functions which have been con- 
structed and applied in recent years mostly in the context 
of dilepton production (with fair success, e.g., in repro- 
ducing the low-mass excess observed by CERES/NA45), 
cf. Ref. [2] for a recent review. This is critical for illu- 
minating the important issue to what extent the relative 
apparent mass shift of about 30 MeV between p-p and 
Au-Au experiments is consistent with such an approach, 
and what responsible mechanisms are. Another question 
that will be addressed concerns total p-meson yields, i.e., 
the p° /ir~ ratio and its evolution from chemical to ther- 
mal freezeout. 

The article is organized as follows. In Sec. II we ap- 
ply microscopic in-medium p spectral functions to com- 
pute thermal 7r + 7r _ emission rates at RHIC, including 
a simple evolution model to assess uncertainties due to 
freezeout conditions. In Sec. Ill we establish a relation 
between thermal emission rates and 7T7t scattering ampli- 
tudes, and employ a realistic 7T7t interaction with (ther- 
mal) in-medium effects to evaluate S-wave (" cr" ) and D- 
wave (/2(1270)) contributions to the emission spectra. 
In Sec. IV we investigate both spectral shape and mag- 
nitude of resonance feeddown to the p-meson yield, as 
well as its evolution from chemical to thermal freezeout 
in central heavy- ion collisions. We also comment on the 
related behavior of the recently measured K* /K ratio 
in the hadronic evolution. In Sec. V we summarize and 
conclude. 
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II. IN-MEDIUM p-MESONS AND tt+tt EMISSION 



B. p — > 7T7T Decays at Freezeout 



A. The p-Meson in Hot Hadronic Matter 

Let us briefly recall the main elements in the calcula- 
tion of the medium-modified p-meson spectral function 
along the lines of our previous work [14-16]. At finite 
temperatures and baryon densities, the p-propagator can 
be cast into the form 
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The in-mcdium sclfcnergy insertions consist of three 
parts: (i) Y, p7T7r encodes the free decay width into 2-pion 
states which in matter is modified by standard "piso- 
bar" excitations, A TV" 1 and NN^ 1 [17-19,14], extended 
to finite total 3-momentum [15] as well as finite tem- 
peratures (most importantly pion Bosefactors), see also 
Refs. [20,21]. Effects of pion excitations into higher res- 
onances, as well as interactions with hyperons, N*- and 
A*-states, are approximated by using an "effective" nu- 
cleon density g e ^ f = (g N + \q b ») [16]; (ii) T, pM de- 
scribes resonant p-interactions with surrounding it, K 
and p mesons (for resonances with established decays into 
p's up to masses ~1.6 GeV) [22,23]; (iii) S p s accounts 
for the resonant p-interactions with surrounding nucle- 
ons, hyperons and baryon resonances [24,14,25,16]. All 
underlying hadronic vertices (characterized by a coupling 
constant and hadronic formfactor cutoff) have been con- 
strained by both hadronic and radiative decay branch- 
ings. In addition, baryonic contributions have been sub- 
jected to a comprehensive fit to photoabsorption data on 
nucleons and nuclei [26]. This, in particular, leaves no 
room for additional t-channel exchanges. 

At collider energies (RHIC and LHC) the net baryon 
densities at midrapidity are rather small. However, as 
pointed out in Ref. [27], CP-invariance of strong inter- 
actions implies that the p meson equally interacts with 
baryons and anti-baryons. Thus, the relevant quantity 
for medium effects is the sum of baryon and antibaryon 
densities, i.e., g B+s = qb+Qb — Qb*(^+P/p), where p/p 
denotes the experimentally measured ratio of antiprotons 
to protons. 

The main medium modification of the resulting p spec- 
tral function, A p = — 2Im£> p , consists of a substantial 
broadening of its resonance shape [16,27], even at mod- 
erate temperatures and, more importantly, baryon den- 
sities. The pole position, on the other hand, is affected 
rather little. This is due to the quite generic feature that 
many-body contributions to the real part of the selfen- 
ergy are of varying sign and thus tend to cancel, whereas 
the imaginary parts are of definite sign and thus strictly 
add up. 



Following Ref. [28] , the thermal rate of two-pion emis- 
sion from neutral p decays per unit four- volume and four- 
momentum can be written as 
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where f p denotes the usual Bosefactor (with chemical po- 
tential n P ), and ImDp depends on temperature, baryon- 
chemical potential p b and possibly other (meson-) chem- 
ical potentials. If the emitted pions reach the detector 
undistorted, the imaginary part of the two-pion selfen- 
ergy is given by the free decay width, — ImE™^(M) = 
m p T v p ac {M), which depends on the invariant p-mass M 
only. Using the standard ptttt interaction vertex, 



£pirir = 9p P 1 " ' (7? X C^T?) , 

its explicit form is (cf., e.g., Refs. [18,14]) 
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with FpTTTr a (weakly momentum-dependent) hadronic 
formfactor (which suppresses the selfenergy at large 
masses M), and k = (M 2 /4 - ml) 1 / 2 the pion decay 3- 
momentum. An invariant-mass distribution is obtained 
upon integrating Eq. (2) over the thermal p 3-momentum 
distribution* , 

x hnD p (M,q;nB,T) . (5) 

If one neglects the (typically weak [16,15]) 3-momentum 
dependence of the in-medium p spectral function, and 
invokes the nonrelativistic limit for the p-meson (scalar) 
density [13], one arrives at (R = dN/d 4 x) 
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The above expression differs from the corresponding one 
in Ref. [13] by (i) an additional 1/M-factor (as well as 
formfactor) in the vacuum width, and (ii) the microscopic 



*if momentum acceptance cuts on the detected pions are 
imposed, the 3-momentum integral in Eq. (5) is subject to 
nontrivial corrections which will then also depend, e.g., on 
the collective flow velocity p-meson at the moment of decay. 
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spectral function, which, most notably, leads to an addi- 
tional M -dependence in lmT, p ot (appearing in the numer- 
ator of the spectral function) , which is reminiscent (but 
not equal) to the free selfenergy. 
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FIG. 1. Comparison of p- meson spectral functions (solid 
lines, in units of GeV -2 ) with invariant-mass distributions 
for p° — > 7r + 7r~ decay rates (dashed lines, arbitrary units) 
in an environment characteristic for p-p (narrow resonance 
curves) and Au-Au (broader resonance curves) collisions. 

It is instructive to illustrate the impact of the vari- 
ous factors in the 7T7t production rate, Eq.(6), relative 
to the p spectral function, for a static environment, i.e., 
at fixed temperature and chemical potentials. To rep- 
resent the p-p case, we use the vacuum p spectral func- 
tion and a temperature of T = 170 MeV characteristic 
for the pertinent hadron production systematics. Com- 
paring the 7T7r invariant-mass distribution with the bare 
spectral function in Fig. 1, one observes a slight shift of 
about -5 MeV, caused by a competition of the thermal 
occupation factor (inducing a downward shift) and the 
various phase space factors (which mostly increase with 
M favoring an upward shift). For the Au-Au case, the 
appropriate conditions correspond to the thermal freeze- 
out stages as extracted, e.g., from blast wave fits to 
hadron p t -spectra. We here adopt the parameters deter- 
mined by hydrodynamic calculations [29] based on suit- 
ably constructed chemical-off-equilibrium hadronic tra- 
jectories [30], yielding T ~ 110 MeV, p^ — 0.5p p ~ 
90 MeV, pn — 365 MeV, etc., translating into a total 
baryon+antibaryon density g B+B = 0.24g - We first 
note that, under these conditions, the p spectral func- 
tion is still appreciably broadened, which, after multi- 
plying with phase space and thermal factors, results in 
a ^-25 MeV relative shift of the maximum in the 7T7r 
invariant-mass distribution, cf. Fig. 1. The rather low 
temperature also induces a significant shape distortion 
which leads to a more pronounced " apparent" resonance 



shift, in the following sense: when evaluating the shift 
of the resonance curves at the full-width half-maximum 
(FWHM), it amounts to about -45 MeV on both sides of 
the resonance. 

For a more realistic description of the emitted (undis- 
torted) pion pairs from the late stages of a heavy-ion 
collision one should allow for a finite emission duration 
over a profile in temperature and density around the ex- 
pected (average) freezeout conditions. A simplified way 
to do this is using a schematic thermal fireball model 
based on cylindrical volume expansion with parameters 
consistent with experiment. The space-time integrated 
emission rate then becomes 

oo 

dN TTr = A ImS^(M) M [ dt V FB (t) P(t,t max ) 

dM 7T° J 

to 

x / ^ n qo ;T(t),p p (t)) 
J qo 

x lmD p (M,q;p B (t),T{t)) Acc(M,q) , (7) 

where Vps(t) denotes the time-dependent (isotropic) vol- 
ume, Acc(M, q) accounts for experimental acceptance, 
and the function P(t, t max ) represents the probability 
distribution of 7r + 7r~ emission around the "average" 
freezeout. For illustration purposes we here assume a 
Gaussian profile, 

P(t, t max ) = P exp [(t - t max f/2al] (8) 

with t m ax the time of maximal emission, a width char- 
acterized by at and a (dimensionless) constant Pq < 1. 
A more detailed treatment with explicit inclusion of pion 
absorption effects will be given elsewhere [31]. In combi- 
nation with the increasing fireball volume, Vp b (t) , the re- 
sulting emission time profile resembles the example given 
in Ref. [11] based on hydrodynamical estimates. The re- 
sulting 7T7T invariant-mass spectra shown in Fig. 2 are 
computed with a t = 1 fm/c and Po = 0.5 (translat- 
ing into an effective emission duration At ~ 1.25 fm/c), 
as well as for two average freezeout temperatures (with 
appropriate pion-chemical potentials [30]) characteristic 
for central Au-Au at full RHIC energy. When employ- 
ing vacuum p spectral functions, the peak positions are 
located at about M„ ~ 760 MeV (i.e., 10 MeV below 
the nominal mass), whereas the in-medium broadening 
effects [16,27] lower the maxima to about M nn — 740- 
745 MeV. Again, the "apparent" mass shift at the center 
of the FWHM is larger due to the distorted line shape. 
One can quantify this feature by defining a " centroid" of 
the mass distribution as 

For the f =110(125) MeV scenario we find m^ nt = 
722(739) MeV and 758(769) MeV with the in-medium 
and free p spectral function, respectively. Furthermore, 
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note that in both cases, the spectral shapes are relatively 
robust with respect to (w.r.t.) variations in the freezeout 
temperature, whereas the yield is reduced by ~20(15)% 
for the free (in-medium) spectral function. We will return 
to the latter issue in Sec. IV B. 
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FIG. 2. Time-integrated tt + tt~ emission spectra from ther- 
mal p°-decays in the late stages of central Au-Au colli- 
sions at RHIC based on Eqs. (7) and (8). The dashed- 
dotted (short-dashed) and full (long-dashed) lines are for 
T = T(tmax) = 110(125) MeV using free and in-medium p 
spectral functions, respectively. 



III. S- AND D-WAVE tttt CORRELATIONS AND 
EMISSION 

A. Microscopic tttt Interactions and Strength 
Distributions 

In addition to direct p-meson decays, a number of 
other sources of (strongly) correlated 7r~7r + emission 
contribute to their invariant-mass distribution. These 
sources are either due to feeddown from resonances (e.g., 
a\ — > pir) or "direct" 7T7t correlations. The former will 
be addressed in Sect. IV; the latter, which we will focus 
on now, are determined by the 7T7r interaction strength 
in a given spin-isospin channel; for an isospin-1 P-wave 
(J 1=11) the strength distribution is essentially saturated 
by the p-meson (i.e., p decays), but non-resonant contri- 
butions are in principle also present. In the S-wave, tttt 
interactions are sizable for both isoscalar (J7=00) and 
isotensor (JJ=02) states. The former represent the well- 
known "cr" -channel, with ongoing controversy on the na- 
ture of the observed resonance structures. For the emis- 
sion strength the relevant quantity is the imaginary part 
of the scattering amplitude, as determined via 7T7t scat- 



tering phase shifts t . In the D-wave, the only significant 
feature is the /2(1270) resonance, with little interaction 
strength below = 1 GeV. 

A microscopic model that accurately reproduces free 
tttt scattering in S-, P- and D- waves up to invariant 
masses of about 1.5 GeV, is the Julich tttt interaction [33] . 
It is based on an effective Lagrangian generating s-, t- 
and it-channel meson exchanges which are iterated to 
all orders by solving an underlying Lippmann-Schwinger 
equation, 



M JI = yjl + yjl M 
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for the scattering amplitude M JI ; V JI and G n7: are the 
two-body potentials (Born amplitudes) and 2-pion prop- 
agator, respectively. In the following we employ a chi- 
rally improved version [34] thereof, which has been ex- 
tended by 7T7r contact interactions to satisfy constraints 
from chiral symmetry (these are necessary for a reliable 
treatment of in-medium effects). In Fig. 3, the solid lines 
represent the imaginary part of the tttt scattering ampli- 
tude in the " cr" , p and /2 channels as obtained from the 
Julich model in vacuum, and including medium effects 
on the two-pion propagator within a hot pion gas [35] 
corresponding to RHIC freezeout conditions (direct res- 
onances, such as pir — ► oi, as well as baryonic effects are 
not accounted for in this Section; a more complete treat- 
ment in " cr" and channels will be reported elsewhere) . 

To obtain the pertinent tt^tt~ emission rates, two fur- 
ther steps are necessary. First, recalling Eq. (2), the 
combination of spectral function and decay width has 



to be related to ImM. Using S 



/.'7T7T 



= v n G 



and 



Mr — VpDpVp (v p : pnn vertex function), one obtains 



ImDp ImSp 



ImA4 n ImG„ 
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with ImG^Tr = k/ (IdmM). Second, the isospin basis has 
to be transformed into particle basis for the two outgoing 
pions. For \II Z ) — 1 10) states (i.e., neutral p- mesons), the 
overlap with 7r + 7r~ pairs is 100%. For S- and Z3-waves, 
however, isospin 1=0 and 1=2 states contribute, 



M: 



r _ = \ M J - I=0 + \m j J= 



(12) 



with only 2/3 of the "cr"- and /2(1270)-mesons decaying 
into tt + tt~ (the other 1/3 into tt°tt°). 

The 3-momentum integrated emission-strength distri- 
butions, dR JI I dM , in the individual spin-isospin chan- 
nels are shown by the dashed lines in Fig. 3 (in arbi- 
trary normalization). Both p- and /2-resonance peaks 



^This is strictly correct only as long as no medium effects 
are included, since the latter may be quite sensitive to the 
underlying microscopic interactions (cf., e.g., Refs. [32] for a 
detailed analysis in cold nuclear matter). 
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are slightly shifted downward mostly due to the ther- 
mal occupation factors, between 10-20(30-40) MeV at 
T=170(110) MeV. The broad structure in the er-channel 
is more strongly distorted by the thermal weights, with 
an appreciable shift of strength towards lower masses, 
especially for the smaller temperature. However, these 
effects are less pronounced than the ones obtained in 
Refs. [11,13]. This is a consequence of the schematic a 
mass distributions used in those works which imply an 
overestimation of the (free) tttt scattering phase shifts, 
and thus of the low invariant-mass strength, in the scalar- 
isoscalar channel. 



B. Thermal 7T ir Emission Spectra 

Implementing the relation (11) into Eq. (5), the spin- 
isospin weighted sum for 7r + 7r~ thermal emission rates 
can be expressed through the tttt scattering amplitudes 
in S-, P- and D- waves as 

dK+iv- k f q 2 dq B 

x ^(2 J + 1) hnMZ+„- (A/; T) . (13) 
J 




0.75 1 

M„ [GeV] 



.5 



FIG. 3. Imaginary part of tttv scattering amplitudes within 
the Jiilich model (solid lines; light: in vacuum, dark: in- 
cluding medium effects at T=110 MeV and ^=90 MeV) 
and corresponding tt + -k~ production rates (arbitrarily nor- 
malized; short-dashed line: T=170 MeV with vacuum am- 
plitudes, long-dashed line: T=110 MeV, ^=90 MeV with 
in-medium amplitudes). 
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FIG. 4. Thermal -k + tv~ emission rates employing the Jiilich 
model scattering amplitudes of Fig. 3. Upper and mid- 
dle panel: using free irn amplitudes at T=170 MeV and 
T=110 MeV, respectively, and ^=0; lower panel: using 
in-medium amplitudes with finite- T pion modifications [35] 
at T=100 MeV and ^=90 MeV. 

The combined results of the Jiilich tttt interaction from 
Fig. 3 are displayed in Fig. 4 for temperatures and pion- 
chemical potentials believed to resemble freezeout condi- 
tions in p-p and central Au-Au collisions. At T=170 MeV 
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(upper panel), where free tttt amplitudes have been em- 
ployed, the spectrum is dominated by the p resonance 
with a slight low-mass shoulder from scalar correlations, 
and a clearly discernible fi resonance. The maximum 
of the p peak is located at M TT = 760 MeV. Reducing 
the emission temperature to T=110 MeV, the stronger 
distortion of the "cr" increases the low- mass shoulder, 
whereas the relative height of the fi peak is suppressed. 
Also note that the p resonance is situated on the decreas- 
ing slope of the "cr" distribution, which entails an addi- 
tional apparent 5-10 MeV downward shift of its position. 
The lower panel in Fig. 4 finally incorporates finite-T 
in-medium effects on the pions in the scattering ampli- 
tudes [35] , as well as finite meson-chemical potentials (as 
inferred from chemistry-conserving thermodynamic tra- 
jectories [30]). One observes a reinforced threshold en- 
hancement, as well as a net shift of the p-peak by about 
-25 MeV, (recall that baryonic medium effects, as well as 
direct p-hadron resonances, are not included here). 

Another noteworthy feature is the minimum structure 
in the Jo (980) resonance region, which is a direct con- 
sequence of the 180°-crossing of the scalar-isoscalar 7r-7r 
phase shifts, i.e., the coherence in formation and decay 
process of the /o (this is different for the " cr" , p and ji 
resonances, which induce the usual 90°-crossing). If, on 
the contrary, a rnr resonance structure is observed in the 
/o(980) region (as, e.g., in pp collisions [9]), it seems to 
require a different production mechanism, such as string 
breaking. More subtle effects associated with the contro- 
versial nature of the /o(980) might also play a role [36], 
e.g., formation through the KK channel. 

C. Final-State Bose-Einstein Correlations 

Apparent mass shifts of p resonances detected via the 
tttt channel are a well-known phenomenon in elementary 
p-p collisions and hadronic jets in e + e~ annihilation at 
high energy. Among the possible explanations [37] are 
interferences between direct and secondary (through tt-tt 
rescattering) p production, as well as phase space dis- 
tortions due to Bose-Einstein correlations between the 
decay-pions and the surrounding ones. Whereas the role 
of the former might be reduced in heavy- ion collisions, 
the opposite can be expected for the latter. 

A naive way of introducing statistical final-state cor- 
relations (which were not addressed in Refs. [10,11,13]) 
into the thermal emission rate consists of modifying the 
free final-state tttt decay by an extra Bosefactor, i.e., 
implementing a factor [1 + 2/ 71 " (tOk',T)] into ImE^^ in 
Eq. (5) . However, this procedure neglects the dependence 
on the total three-momentum q of the pion-pair (w.r.t. 
the thermal frame). A better approximation is obtained 
by performing an angular average according to [38] 

(l + r(k;T) + r(-k;T)) = 



+i 

= i J d cos6 [l + r(k' + ;T) + r(k'_;T)} 
-l 

T l f sinh(^[7K+/3fc)-^]) ) 
0jk \sinh(^[ 7 K-/3fc)-^])/ 

which depends on the pair 3-momentum and invariant- 
mass through the Lorentz factors (3—q/qo and r )=qT)jM^ 
(k'± are the pion momenta in the thermal rest frame). 
Consequently, Im£ p7r7r in Eq. (5) has to be evaluated in- 
side the momentum integral. An upper estimate of the 
final-state Bose correlations may be obtained by insert- 
ing the same temperature and chemical potentials as in 
the spectral function in Eq. (5). This leads to an appre- 
ciable impact on the tt~tt~~ emission spectra, cf. Fig. 5: 
in the p-resonance region, the rate increases by 20 % 
("stimulated" emission), accompanied by a small down- 
ward move (5-10 MeV) of the peak position. Towards 
the two-pion threshold, a more pronounced " cr" shoulder 
emerges. The finitc-g corrections moderate the Bose en- 
hancement close to the two-pion threshold, most notably 
for low temperature and large pion-chemical potential. 
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FIG. 5. Thermal tt tt~ emission rates employing the Jiilich 
model scattering amplitudes with additional account for fi- 
nal-state Bose-Einstein enhancement. Dashed-dotted lines 
represent the results of Sec. IIIB (solid lines in Fig. 4), 
whereas solid and dashed lines include the final-state Bose- 
factors with and without finite-g corrections, respectively. 
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IV. RESONANCE FEEDDOWN AND RATIOS 

We now return to the question of feeddown contribu- 
tions to 7r + 7r~ distributions, concerning both their mag- 
nitude and spectral shape. For simplicity, we will focus 
on the p meson, which is well-established in a large num- 
ber of mesonic and baryonic resonance-decay branchings 
(less so for the scalar "<r"- and tensor / 2 -mesons). 



little impact either. Final-state Bose enhancement fac- 
tors, implemented along the lines of Sec. Ill C, can lead to 
net peak-shifts of maximally -15 MeV. Note that the rel- 
ative contribution from resonance feeddown is expected 
to decrease with decreasing temperatures. 



B. p/ty and K* / K Ratios 



A. Spectral ir-iv Shape from ai Decays 

As a typical example for resonance feeddown, we here 
concentrate the ai(1260) — ► pir decay (which gives the 
largest individual contribution). Starting point is the 3- 
momentum integrated emission rate, Eq. (6), which for 
the present case takes the form 



dRai — >piv 6 

dM a ~ 7T 



Im£ aip7r (M a ) lmD ai (M a ;p B ,T) 



\^f) 3/2 e-^-^' T . (15) 



The aipir selfenergy is related to the pertinent partial 
decay width via — ImS aip7r = M a T ai ^ p7T . Including the 
p-mass distribution through its spectral function A p (M), 
one has 



M„ 



i r MdM 
7m 2 J ~ 

2m, 



12ttM2 



A P {M) q c 



(16) 



where q cm denotes the three-momentum of the decay 
products p (mass M) and it (mass m w ) in the a\ rest 
frame, v p - !Tai the vertex function, and M max = M a - m n . 
Using the differential form of Eq. (16), and integrating 
Eq. (15) over the ai-mass distribution, yields the p-likc 
7r + 7r~ production rate from a\ deacys as 



dR, 



"'-*(* + *-) =_£ / dMn M n 



d 3 k M a dT ai ^ p 



dM tt J " U J (2tt) 3 k dM 

x f ai (k ;n ai ,T) lmD ai (M a ,k;p B ,T) , 

(17) 



where we reintroduced the explicit integration over the 
ai Bose function, f ai , with fc = (M 2 + k 2 ) 1 / 2 . An im- 
portant ingredient in evaluating Eq. (17) is the a\ spec- 
tral function; we use a microscopic description with a pir 
selfenergy which gives an accurate description of recent 
axialvector r-decay data. 

It turns out that the resulting 7r + 7r~ mass distributions 
are not much affected (downward peak-shifts of 5-10 MeV 
even at low temperatures), which is primarily due to the 
fact that the free p spectral function enters into the ex- 
pression for the ai-decay width, Eq. (16). Broadening 
effects in the underlying a\ spectral function have very 



For a thermal hadron gas at T = 180 MeV in full 
chemical equilibrium, an additional <~ 60% of the equi- 
librium p-number will emerge from resonance feeddown 
(this excludes subthreshold contributions of the type 
w — > pn or A r (1520) — > Np which do not feed into 
the p- meson peak); the largest component of about 13% 
each arises from ai(1260) — > pn and a 2 (1320) — > pn de- 
cays (at p,B = 0, about 10% arise from anti-/baryon de- 
cays). The resulting p°/ir~ ratio ranges between 0.106 
and 0.125 depending on whether weak decays arc in- 
cluded in the tt~ yield or not. Such values are in good 
agreement with previous measurements in e+e~ and pp 
collisions at ^/s ~ 10-100 GeV. The recent (prelimi- 
nary) STAR measurement for y/s = 200 GeV pp reports 
p°/ir- = 0.173 ± 0.026 [9]. This value, however, does 
not account for the possibility that other n + n~ sources 
(most notably "<x" correlations) may "contaminate" the 
p resonance region. From the upper panel of Fig. 4 one 
infers that this could lead to a up to 20% correction of 
the p° yield. 

An important question is how the p° /ir~ ratio devel- 
ops when going to (central) heavy-ion collisions. The p 
meson constitutes one of the strongest agents for pion in- 
teractions, an thus for maintaining thermal equilibrium. 
This means that p <-» nn stays in relative chemical equi- 
librium, and therefore p p = 2[i n . Finite p^s build up 
between chemical and thermal frcczeout due to effective 
pion-number conservation (and likewise for kaons, etas, 
and even antibaryons), cf. Ref. [30] for a hadronic evolu- 
tion (i.e., thermodynamic trajectory) under RHIC con- 
ditions. To obtain the total number of p° and ir~ in the 
later stages one again has to include the feeddown cor- 
rections of all the resonances, but with their respective 
(effective) chemical potentials (e.g., p ai = p p +p^ = 3p, v , 
etc.). Starting the evolution from chemical freezeout 
at ((j, n ,T) = (0,180) MeV, the p°/7r" ratio is first al- 
most stable, being reduced by only 15% at (//„., T) = 
(60, 140) MeV. It then starts falling faster being reduced 
by 30% at (^,T) = (85, 120) MeV, i.e., p°/-K~ = 0.075- 
0.087, and to 0.05-0.06 at (/i w ,T) = (100,100) MeV. 
These ranges thus characterize (equilibrium-) conditions 
close to thermal decoupling in central Au-Au collisions. 
The (preliminary) STAR measurement in peripheral Au- 
Au lies significantly higher, at 0.183±0.027 [9]. Again, 
7T7T S*-wave contributions in the spectra could lower this 
value. Also, the actually detected number of p° — > ir + ir~ 
decays might be affected by an interplay between (pion-) 
absorption and (p -) regeneration in the vicinity of ther- 
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mal freezeout. Here, transport simulations should be 
an adequate tool, which find a p°/tt~ ratio of around 
8% [39]. 

Let us briefly comment on a closely related quan- 
tity that has recently been measured at RHIC ener- 
gies, namely the K*(892)/K ratio. For 200 GcV p-p 
collisions STAR obtained 0.386±0.029 (preliminary), to 
be compared to a statistical model value of 0.35-0.36 
at T=180 MeV. For 0-10% (0-14%) central Au-Au at 
200 (130) AGeV, the ratio was found to decrease to 
0.205±0.033 (0.26±0.07) [40,41]. On the one hand, it has 
been argued that this implies, based on the vacuum K* 
lifetime of about 4 fm/c, an upper limit on the duration 
of the hadronic phase, which follows if one assumes no re- 
generation. If, on the other hand, the reaction K* «-> Kir 
is assumed to maintain (relative) chemical equilibrium, 
one has uk* = P-k + Pk- Along the thermodynamic tra- 
jectory discussed above one then finds the K* / K ratio 
to slowly decrease, passing through values of 0.235-0.295 
at temperatures T=100-120 MeV (and ^ mes >0). Thus, 
even complete regeneration (which, in the denser stages, 
could be facilitated by in-medium broadening effects) 
cannot be ruled out by the current measurements. The 
same pattern emerges for the K°* jh~ ratio. The prelim- 
inary STAR value for 0-20% (0-14%) central Au-Au at 
200 (130) AGeV is 0.033±0.004 (0.042±0.011) [9,41]. For 
the hadron gas value at chemical freezeout (T=180 MeV) 
we find K°*/h~ =0.041, in good agreement with the 
statistical model values of Ref. [42]. Around thermal 
freezeout (T=100-120 MeV) the equilibrium ratio has de- 
creased to K°*/h- =0.028-0.035. 

V. SUMMARY AND CONCLUSIONS 

In the present article we have investigated spectral 
shapes and yields of 7r + 7r~ invariant-mass spectra aris- 
ing from thermal sources characteristic for the decoupling 
stages of high-energy heavy- ion collisions. Our analysis is 
based on realistic tttt interactions in S-, P- and D- waves, 
adjusted to free scattering phase shifts. 

For the P-wave, in line with previous studies [11,13], 
thermal phase-space factors were found to induce up to 
10 MeV downward shifts of the p-peak position when as- 
suming the vacuum resonance profile. However, when us- 
ing microscopic in- medium p spectral functions, for which 
a broadening of about 80-100 MeV in dilute hadronic 
matter has been predicted earlier, the combination with 
thermal 7T7t phase space leads to a peak-mass reduction 
of about 30 MeV. In fact, low-mass tails of the tttt distri- 
bution lower its centroid by up to another 20 MeV. These 
findings arc more pronounced than the corresponding 
ones of Ref. [13] (based on schematic spectral functions). 
Although the effects of in-medium modifications qualita- 
tively reflect the relative changes between p-p and Au-Au 
collisions as extracted from preliminary STAR data [9], 
the experimentally observed absolute peak shifts of about 



-40 MeV and -70 MeV, respectively, are not reproduced 
by our results. 

For S'-wave emission, especially in the scalar-isoscalar 
("er") channel, we emphasized the importance of employ- 
ing realistic tttt scattering amplitudes (strength distribu- 
tions), which sensitively respond to distortions by ther- 
mal occupation factors close to freezeout. This is manda- 
tory for a reliable assessment of in-medium effects on the 
"ct", which, as for the p, can be related to partial chiral 
symmetry restoration. We found that thermal medium 
modifications and pion Bose-enhancement lead to an ex- 
tended threshold maximum (M^^ =300-500 MeV) in the 
S'-wave 7r + 7r~ emission rates. However, its magnitude 
relative to the P-wave contribution is moderate, trans- 
lating into a low-mass shoulder of the p resonance. In ad- 
dition, the decreasing " a" slope in the p- resonance region 
induces another -10 MeV apparent shift of the observed 
p peak. The general trend in the preliminary STAR data 
which indicates a decrease in the height of the p-peak 
relative to the low-mass continuum when going from p-p 
to Au-Au collisions, is qualitatively reproduced by our 
results. To what extent this reflects a p broadening, and 
what additional mass shift is required, should be evalu- 
ated by a more quantitative comparison to experimental 
spectra. 

Finally, we have assessed total p° yields, includ- 
ing feeddown corrections. We found that the hadro- 
chemical freezeout value for the p°/tt~ ratio (^11% at 
T=180 MeV) decreases by only about 30-40% towards 
thermal freezeout (T=110-120 MeV), due to large pion- 
chemical potentials. 

For future investigations, a more complete treatment 
of in-medium effects on the " cr" , combining finite temper- 
atures and densities, is desirable. The /o(980) — > tt + tt~ 
contribution clearly deserves further study. For a realistic 
comparison with experimental data, Dalitz decays (e.g., 
ijj — > (7T + 7r~)7r°) have to be included. Also, pion absorp- 
tion effects need to be addressed explicitly In addition, 
the relation to resonance modifications in high-energy pp 
and e + e~ collisions ought to be understood better [37]. 
Only an overall consistent treatment will eventually al- 
low for a reliable extraction of in-medium effects from 
short-lived resonance spectroscopy in the heavy-ion en- 
vironment. 
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